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Cosmological Models 
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Abstract We study two cosmological models: A non- 
minimally coupled scalar field on brane world model 
and a minimally coupled scalar field on Lorentz invari- 
ance violation model. We compare some cosmological 
results in these scenarios. Also, we consider some types 
of Rip singularity solution in both models. 

Keywords Scaler- Vector- Tensor Theories, Theories of 
Extra Dimensions, Brane world Cosmology, Late-time 
Acceleration, DGP Brane Cosmology, Lorentz Invari- 
ance Violation. 



1 Introduction 

In theories of extra spatial dimensions, ordinary mat- 
ter is captured on the brane but gravitation extends 
through the entire space-time (lAr kani-Ha med et al. 
1998t iRandall k Sundrumlll999HDvali et al.ll2000l) ~In 
these scenarios, the cosmological evolution on the 
brane is taken by an effective Friedmann equation 
that combines the effects of the bulk in a n on-trivial 
kind (|Binetruv et alJfeoOOallbl ; lMaartensll2004l ). In other 
view point, the important result in brane models is 
an alternative scenario for late-time expansion of the 
universe. This result predicts deviations from the 4- 
dimcnsional gravity at limited distances. In other 
hand, the model considered by Dvali, Gabadadze 
and Porrati (DGP) is different from above models 
since it also predicts deviations from t he standard 4- 
dimen s ional gravity in large di stances (jBinetruv et al 
2000at ICollins k Holdoml |2000h . Generally one can 
study the effect of a caused gravity term as a quantum 
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discipline in any brane world model. The existence of a 
higher dimensional embedding space lets for the bulk or 
brane matter that can evidently affect the cosmological 
evolu tion on the brane( L~anglois k Rodriguez-Martinez! 
2001J). A special form of bulk or brane context is a 
scalar field. Scalar fields take an important key both in 
models of the early universe and late-time acceleration. 
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2003 



In other hand, Lorentz invariancc violation models 
(LIV) h as been considered in the scalar-vector-tensor 
theories ( Kanno k Soda 2006f ). It has described that 
Lorentz violating vector fields influence the dynamics 
equations in the inflationary models. An interesting re- 
sult of this model is that the exact Lorentz violating in- 
flationary solutions are depended on the absence of the 
inflation potential. In this case, the infla tion is exactly 



collaborated with the Lorentz violation! Arianto et al 



20071 ). Therefore, we study this symmetry breaking on 
the dynamics of equation of state for some cosmological 
aspects. 



1999; 



Riessetal 



Some evidences from su pernova data (IPerlmutter et al 
" 19981). CCMB) results dMiller et~ 



19991 Ide Bernardis et alj|2000t lllanany et al.ll2000l) and 



(WMAP) data (ISpergel et all 2005 ; iPage et al 



2003 



Spergel et al]|2007t iReiss et alJbooi IZhao et alj|2012h 

point out an accelerating phase of cosmological expan- 
sion and this characterize shows that the picture of 
universe by pressureless fluid is not enough; the uni- 
verse should contain some type of additional negative 



pressure known as dark ener 



in this field see (Arianto et al 



;y(For brief introduction 



20071 )). Also, the merged 



analysis of the WMAP data with th e supernova Legacy 
survey (SNLS) (|Spergel et al.1 120071 ). compels the equa- 
tion of state Wde, in accord with 74% donation of dark 
energy in the currently accelerating universe. More- 
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over, observations show so me sort of a dar k energy 
equation of state, Wde < — 1 (jReiss et al.ll2004l ). Hence, 
a practical cosmological model should accept a dynam- 
ical equation of state that may have crossed the value 
Wde = — 1, in late time of cosmological evolution. 
In following, we consider cosmological results of a non- 
minimally coupled scalar field on the brane and a min- 
imally coupled scaler field in LIV model. Also we de- 
termine late-time behavior of our equations and obtain 
some restriction on the parameters of models to have 
an accelerating universe. 

We understand that dark energy is an important prob- 
lem in modern cosmology, especially, if the equation 
of state parameter u less than — 1. In this article, we 
study basically solving Friedman equations and obtain 
the evolution of the effective equation of state param- 
eter uj. The phantom crossing conditions that we ob- 
tained for these set-ups, themselves are interesting at 
some degree. Furthermore, we know that just consid- 
ering background quantity is not enough in such kind 
of unconventional cosmological models. The current 
constraint on u) is obtained by combining the result 
of the observation of CMB, which means that the in- 
formation of the evolution of perturbations should be 
also included. Also we understand in the brane world 
model, the perturbations of brane are coupled with 
these of bulk, which gives nontrivial effect. Therefore, 
in first stage we just study the background dynamics 
and the evolution of perturbations will study in future. 
In other view, there is a no-go theorem which explic- 
itly points out a conventional dark energy model in- 
volving single degree of freedom in the frame of stan- 
dard Einstein gravity is forbidden to realize such a 
scenario. This no-go theo rem was proven in the ap- 



pendix of (jXia et al.l 120081 ) . In this regards , we point 



out the original work in ( Feng et al. 20051) and sev 
eral papers addres s ing this scenario (jCai et al] 12007 ; 
Cai fc Wan j l2008t ICai et al.l |2010|) . This theorem 
is exactly the reason why we study a number of 
nonconventional dark energy mode ls in realizing the 
Quintom scenario ()Zhao et al.l 120121 ) . such as the non- 
minimally coupled on brane world model and the 
Lorentz-violating model considered in the present pa- 
per. However, a non-minimally coupling scalar field 
identified on the brane model and scalar field coupling 
minimally to gravity in LIV model have some simi- 
lar cosmological results, this is an interesting theoret- 
ical motivation of these models. Also we study other 
solutions admitted as Rip singularity, that occur in 
the condition uj < — 1 increases rapidly. However, it 
possible different types of singularity, depending en- 
ergy density and scale facto r how increases with time 
( Frampton et al.ll201ll . |2012|) . 



2 Non-minimally coupled scalar Field on the 
Brane 

Here we study a brane world model where a scalar 
field is coupling non-minimally to the Ricci scalar of 
the brane. In following we only consider a scalar field 
in the matter Lagrangian without taking into account 
baryons, cold dark matter, and radiation. This kind of 
analysis suffers from a potential maybe cause that the 
model not be able to explain the evolution of a realistic 
universe at background level when confronting with ob- 
servations. However, this problem can solve by suitable 
fine-tuning parameters of model. 

The action in the absence of ordinary ma t ter can 
be g iven as( Bouhamdi-Lopez fc Wands 2005 ; Mvunj 
200ll) 



S 



(1) 



where we have made a general non-minimal coupling 
a(<f>). For simplicity, in following we set k^ = 8ttGn = 
1. We can obtain Einstein equations with variation of 
the action respect to brane metric 



(2) 



where Tjj,u , energy- momentum tensor of the scalar field 
non-minimally coupled to gravity is taken by 
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Ua(4>)- V M V i ,a(</>), 



(3) 



where □ shows 4-dimensional d'Alembertian. We study 
the FRW universe with line element as following 



ds 2 = -dt 2 + a 2 (t)dS fc 5 



(4) 



where gE/. is the line element for a constant curva- 
ture k = +1,0,-1. The Ricci scalar obtain with the 
equation of motion for scalar field as 



V^^ = V -a'R[g], 



(5) 



where a prime denotes the derivative of each parameter 
with respect to 0, that may be written by 



a ■ dV , „. , 
■3-4>+-rr=a'R\g]. 
a dip 



(6) 



where a dot means the derivative of each parameter by 
respect to t The intrinsic Ricci scalar for a FRW brane 
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give as 

R[g] = 6^H + 2H 2 + ^J, 

and Fricdmann's equations arc determined by 

a 2 k p 
a 2 a 2 3 ' 

and 



a 1 . 

- = — e (p + 3p). 
a o 



(7) 



(8) 



(9) 



We take a scalar field, <j>, only depended on time. So, 
with Eq(3), we obtain 



p = a- 1 ( ^ 2 + V{4>) - 6a'H(j) 



(10) 



p = a- 1 \^-<p 2 -V(<j)) + 2(a't + 2Ha'<p + a"<p 2 )j (11) 

where H = - is Hubble parameter. Now equation of 
state has the following form 



w 



P 
9 



ft - 2V((f>) + i(a'(f> + 2Ha'4> + a"ft^j 
ft + 2V{<t>) - 12a'H<j) 



(12) 



When <j) = 0, we have p = —p. In this illustration p 
depended on a and V((f>), it has the duty of a cosmo- 
logical constant. In the minimal case which ft < V (</)), 
by Eq(9), we take p < — | that shows a late-time ac- 
celerating universe. In non-minimal case, the position 
depends on the choice of non-minimal coupling param- 
eter. In following we show how for a suitable range of 
coupling parameter, a late-time accelerating expansion 
can be described. We first consider a moving domain 
wall of brane world to discuss quintessence behavior, 
then we study a special non-minimal coupling for late- 
time acceleration. 

2.1 Late-Time Acceleration in a Brane world Model 



conditions ([Carroll 2004). Here we study the following 
line element (|MvunsJ 120011 ) 



ds 



5± 



k 



rj± 



( 2 



dt 2 



r 2 ~ i'i 



-dr 2 



r *fijdx l dx 3 , 



(13) 



where ± is for left(— ) and right (+) side of the moving 
domain wall, also £ is curvature radius of AdSs mani- 
fold and "fij is the horizon metric. rj± =/= cre ates the 
electric part of the Weyl tensor on two sides! Mvunj 



200lh . 

We assume usual mater on the brane has a perfect fluid 
form, T M „ = (p + p^u^Uy + ph^ where p = p m + a 
and p = p m — a . Energy density of the confined mat- 
ter on the brane and analogous pressure are shown 
with p m and p m , while a is the brane tension, with 



Israel junction conditions (jCarrollI 120041 ) and Gauss- 
Codazzi e quations, we obtain generalized Friedmann 
equations (|Mvungll2001l ) as following 



k 



V 



3G 



(14) 



'- = -^(l + 'Sw)-^-^p 2 n (2 + 3w) (15) 
a 6 cp 3d 

where we require a ^-symmetry by rj + = r\- = rj and 
we have supposed p m = wp m . following we take case: 
rj = 0. For 77 = 0, each sub-manifolds of bulk space- 
time are accurate AdSs space-times. Now we study 
a confined non-minimally coupled scalar field on the 
brane and consider cosmological aspects. We use energy 
density and pressure of scalar field (10) and (11) and 
by equation (15), hence for cosmic expansion have 



a 
a 



1 

Go: 



2 + V{4>) - 6a'H<j) x 



^ 3 ^ 2 -2V{4>)+A(a'4> + 2Ha l 4) + a"4> 2 ) 
<p 2 + 2V(<f>) - l2a'H<j) 



We study a bulk config by two 5-dimensional anti de 
Sittcr-Schwarzschild (AdSSs) black hole spaces com- 
bined on a moving domain wall. To insert this mov- 
ing domain wall into 5-dimensional bulk, it is required 
to indicate normal and tangent to the domain wall 
by determination of normal instructing to the brane. 
We take that domain wall is identified at coordinate 
r = a(r) where a(r) is considered by Israel junction 



e 2 (1 ; 



-(f> 2 + V{4>) - 6a'H<j) x 



36a 2 V 2 



3 2 - 2V(4>) + 4 (qj + 2Ha'<p + a" ft) ' ) 



(f) 2 + 2V{(j>) - 12a'H<p 

where H — - is Hubble parameter. This is a complex 
equation, to describe cosmological aspects, we must 



•1 



study some limiting cases or specify a{4>), V(4>) and 



Equation (16) shows, the condition for an accelerating 
universe depended on the choice of non-minimal cou- 
pling and the scalar field potential. Hence we determine 
the following non-minimal coupling 



(17) 



Fo r an accelerating universe we shou l d have Pm+3p m < 
|Bouhamdi-Lopez k Wands! 120051 iFaraonilboOOl) . If 
we take V > 0, we obtain 

(1 - 30«^ 2 - V{4>) - 3£t/>(4> + H4>) < (18) 
dynamics of this scalar field is given as 

4> + 3H<j) + ^R(f>+ — =0, (19) 
d(p 

we obtain 

(l-30^ 2 -F(^) + 3£ 2 /ty 2 +6£W+3£^ < 0. (20) 

d<p 

With p m in this case, we determine 

y=(l-^ 2 )p m -2V( ( f>) 
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-3£ 



3£ 2 i?0 2 



3^<0 



(21) 



that is the requirement to have an accelerated universe. 
In following, we study the weak energy condition p m > 
and limit our consideration to the case by £ < 1/6. 
Therefore we obtain 



-2V + 3£</>^ < y < 



and a main condition for cosmic acceleration is 

1 



6 



(22) 



(23) 



For the model have cosmic acceleration by £ > 0, the 
potential V(cj>) must change with <fr slower than power- 

2 

law potential V c (<f>) = Vo\4^j ^ (Fig !)• In other hand, 
when £ < 0, the main condition for cosmic acc eleration 
needs V grow faster than V c as <f> increases (jFaraoni 
2000). Now we take a simple special example to show 
how this model process. By potential of the form 
V{(j)) = condition (23) takes A^l - ^) > 0. 

In this case just for £ < 2/3n we have the accelerated 
expansion. 




Fig. 1 V vary for different values of the <j> and Non-minimal 

2 

coupling parameter £ for V = (f>^ . 

More general case is r\ =^ in equations (14) and 
(15). In following we convert the generalized Friedmann 
equation as 

X + 4- \x + t^- (1(1 - 2 O0 2 - V{4>) - 2^U + 3-< 
a 12a V 2 V a 



= 0, (24) 

where \ = is dark energy field. Then we take 

4>(t) w where A is a constant and /3 an unspecified 
power will be studied. For scale factor we take the 
ansatz a(t) sa Bt v where B is a constant and v to be 
determined. For scalar field potential we take V(<p) = 
A^-. Therefore we obtain 



a 
a 



7T- + ITT—Pi 1 - 3u>) + 7T- (V(0) + A) 

12aw 24au^ ; 6a V / 



12a 



and 



(1 + 4a") 



2u 



a 
3a" 



(25) 



■ + io^(n0) 

u u \ 



A 



+4 a'(u + aV)^ (i _ + 8a> + aV) ff2 (2g) 



where u = 6 



i-ea-fo^- 



We take usual matter 



on the brane has an equation of state such as p = 



5 



Hence we have 




t 2P 



a(l-3w)D\ 1 
6 J 



4^D(l-3tt)(|-fV. 



1 



2/3 



(l-80+4£(3i/-i; 



and scalar field equation given as 
A XA 



(27) 



t P+2 



3 ? t 2 / 3 



5£ A 

3 t* 3 



A 



A i, 16 .\ A 2 

2 + |( 1 "T f ) A M 



(28) 



where /3 > 1. Two conditions must solve numerically 
to limit the values of v and /? for an accelerating ex- 
pansion. After some calculations we obtain 



(1 - Zw)ja 
9 



±1 
12 



D = A/3(p + 1 - 3u), 



(29) 



(30) 



where we passed over terms of order 0(t 3,3 ) and 
higher. Now, if we arrange A, since v(2v — 1) = 



<t(1 — 3w) t-a 

— i we find 



144 



where 



(3^-l) 2 



16£i/(2i/ - 1) 



(31) 



^jq-^y. For v < 1 the parameter /? is real 
for any value of £ , for v > 1 the exponent /3 is real only 
if we have 



1 (3*/-l) 2 
16 2^(2^- 1)' 



(32) 



These equations limited the values of non- minimal cou- 
pling £ for obtain an accelerating expansion(Fig 2). 
Hence, it seems that non-minimally coupled scalar field 
identified on the brane world supplies natural candidate 
for late-time expansion. 



1.8 2.0 




Fig. 2 Variation of £ respect to v with equation (32) (up) 
and Variation of ft respect to v and £ with equation (31) by 
sign +, for sing — in this equation (31) /3 has negative value. 
We choose special range of parameter £ with equation (32) 
in (down) figure. 



Also we can use equation (12) for determine dynamical 
Equation of state to(t). 

Figure (3) shows possible crossing of phantom di- 
vided barrier in this framework, this result is important 
because previous considerations have shown, crossing 
phantom divided line with a scalar field non-minimally 
coupling gravi ty exist (jNesseris fc Perivolaropoulosi2007t | 
Vikmanl[2005h . 



In following, we explain the physical reasons why the 
crossing of the phantom divide can occur in this model. 
Here, we have a scalar field that non minimally coupled 
in gravity and an ordinary matter. Ge nerally, according 
to ( Nesseris fc Perivolaropoulosl 2007 ). these sources of 
energy-momentum can switch together and describe 
evolution of universe. In this regard, an encouraging 
candidate for the dark matter is either a positive cosmo- 
logical constant or a slowly developing scalar known as 
quintessence. Also the quintessence can be determined 
as an alternative method to solve the cosmological con- 
stant problem. This is possible because instead of the 
fine-tuning, the quintessence gives a model of slowly 
decaying cosmological constant. However, there exists 
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0.06 
0.04 xi 



Fig. 3 Dynamical Equation of state aj(t) vary for different 
values of the t and Non-minimal coupling parameter £ by 
equation (12) and specific choose for parameters space. 



another question for the quintessence now. This is to 
inquire whether the expansion will keep accelerating 
forever or it will slow down again after some time. This 
is very similar to the leaving problem in the inflation- 
ary cosmology. According to the theory of quintessence, 
the dark energy of the universe is controlled by the po- 
tential of a scalar field <f> which is still curling to its 
minimum at V = 0. 



2.2 Rip Singularity 

The new interest to the models with the crossing phan- 
tom divi de w < — 1 is the ir prediction of a Rip sin- 
gularity ( Starobinskv 2000t ). Generally, the scale fac- 
tor of the universe gets infinite at a finite time in the 
future which was named Big Rip singularity. There 
were proposed some scenarios to improved the Big 
Rip singularity: (I) To determined phantom acceler- 
ation as temporary phenomenon. This is a number 
of scalar potentials. (II) To explain for quantum ef- 
fects which maybe de lay/stop the singularity happening 
|Elizalde et alj|2004h . (Ill) To change the gravitation 



itself in such a way that it seems to be observationally- 
fricndly from a side but it restores to singularity. (IV) 
To couple dark energy with dark matter in the spe- 
cific way or t o use specific form for dark energy equa- 
tion of state (jBamba et al.l 120081 ). Pay attention to for 
quintessence dark energy, other finite-time singularities 
may appear. For instance, type II (sudden) singularity 
or type III singularity appears with finite scale factor 
but infinite energy and/or pressure. The near exami- 
nation shows that the situation u < — 1 is not enough 



for a singularity event ( No iiri et al. 20051) . First of all, 
a transient phantom cosmology is possible. However, 
one can make such models that w asymptotically tends 
to -1 and the energy density increases with time or 
corpses constant b ut there is no finite-time singularity 
(jNoiiri et al.l 12005). Of course, most apparent case is 
when Hubble rate inclined to constant (cosmological 
constant or asymptotically de Sitter space), which ca n 



also point out the pseudo-rip ( Frampton et al.l 12012 ). 
most intere sting situation is conne cted with Little Rip 
cosmology (IFrampton et al.l 120111 ) where Hubble rate 
go to infinity in the infinite future. The goal point is 
that if lu approaches -1 adequately rapidly, then it is 
possible to have a model in which the time needed for 
singularity is infinite. However, this situation does Not 
occur in our model, because H is finite. 
In following, if we want to determind some types of 
Rip singularity, we must use Hubble parameter, energy 
density and pressure equations that obtained in above 
equations. Therefore, after calculations with £ = 0.14 
(figure 4), we take results as following: 

case 1 : type II singularity (sudden) 
conditions: scale factor and energy density arrives finite 
value but p —t oo 

In this model, we find that scale factor, pressure and 
energy density are finite. Therefore, the model does not 
predict a sudden singularity. 

case 2 : type III singularity 
Conditions, scale factor arrives finite value but (p — > oo) 
and | p |— > oo 

This means the energy density becomes so rapidly 
with time and scale factor does not arrive the infinite 
value. The goal difference between case (1) and (2) 
for this model is that the potential of scalar field has 
a stake in a case of singularity. Here, according nu- 
merical calculations, if scale factor arrives finite value, 
energy density and pressure are finite. Hence, type III 
singularity does not occur in this setup. 

case 3 : pseudo-rip singularity 
conditions : Pressure p — > —p in infinite time and cause 
to H -> 

This means the expansion of the universe asymptoti- 
cally appeals to the exponential area. Obviously in our 
model with above setup, pseudo-rip singularity predict, 
because in infinite time, when H arrive to zero, p — > —p. 



3 A Lorentz Violating Model 



With using dKanno fc Sodal bood lArianto et ail 12007 : 
Sadatianl 120121 ). we consider the cosmological aspects 
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\V^u v V,u v - /3 2 V W.u, 



-fa (V^f " fau* i u v V^u a V v u a + A + 1)] (35) 



S,i. 



d A x, 



-9 £-4 



(36) 



This action has a non-gravitational degrees of free- 
dom in the structure of Lorentz violating scalar-tensor- 
vector theory. We take u^u^ — —1 and the ex- 
pec ted value of vector field is < OIm^m^O > = 
-1 dKosteleckv fc Samuell 1 1986ft . (i = 1,2,3,4) 

are unrestricted parameters by dimension of mass 
squared , also is th e Lagrangian density for the 
scalar field in this model (|Kanno fc Sodall2006l ). If uni- 



verse be homogeneous and isotropic, we consider the 
universe by metric as 



ds 1 



-N 2 {t)dt 2 



e 2ait) S ij dx i dx j 



(37) 



where Af is a lapse function and the universe scale i s 
given by a jKanno fc Sodalbood lArianto et aTll2007ft . 
If the action vary with respect to metric and selecting 
a fitting gauge, field equations obtain as 



Fig. 4 Variation of pressure(up) and energy density(down) 
respect to time for a non minimal coupling model(where 
£ = 0.14) 

of Lorentz invariance violating model. We empha- 
size this model parameters are tightly bounded by 
both astronomic al and cosmological experiments, some 
works such as ( Arianto et al. 2007 . 20081 ; Zen et al 
20091: lArianto et al.ll201oF explained, the model param- 
eters are able to be bounded by comparing the model 
with cosmological observations. Also a global anal- 
ysis and a group of canonical values for the model 
parameters to constrain th ese parameters given in 
(jKosteleckv fc Samuellll986l ). so that the model is ex- 
pected to be viable and background evolution is rea- 
sonable. 

We want to obtain a connection between Lorentz Invari- 
ance violation parameter and dynamics of scalar field. 
First we define an action for a representative scalar- 
vector-tensor theory which allows Lorentz invariance 
violation as 



S — Sn "T" S u -\- Ss , 



(33) 



where the actions for the vector field S u , the tensor field 



S g , and the scalar field Stj, are assume as 



So 



d 4 Xy 



1 

16ttG 



R 



(34) 



(38) 



where T, 



rp{u) 



T$ is the total energy-momentum 
tensor, and energy-momentum tensors of vector and 
scalar fields given by T^t . The time and space ele- 
ments of the total e nergy-momentum tensor obtain as 
( Arianto et~ai]|2007ft 



71U 
-L n 



Pa + Pi, 



(39) 



where the energy density and pressure of the vector field 
determine as 



p u = -3PH 2 , 

Pv 

P = fa + 3 fa + fa 



H' fa\ 
3 + 2— + 2^ ) faH 2 , 



(40) 
(41) 
(42) 



a prime point to the derivative of any parameters by 
respect to a and H = da/dt = a is the Hubble 
parameter. The energy equations for scalar field, <f> and 
the vector field u given as 



p' u + 3{Pu+Pu) = +3H 2 /3' 



P'^ + HP4,+P4>) = -3H 2 P' , 



(43) 



(44) 
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So, the total energy equation for two the vector and the 
scalar fields write as 



p' + 3(p + p) = , = p„ + , 



(45) 



Now, dynamics of the model determine with the 
following; Friedman n equations (|Kanno fc Sodal 12006 ; 
Arianto et aj]|2007h 



1 



1 



8ttG/3 
1 



H 2 - 1 , 



(46) 



+ ^)-H 2 ^.(47) 



8irGf3 J 



1 

6^ p 



For the scalar section of this model we assume the fol 
lowing Lagrangian 

'/ 



C 



(48) 



where (V^) 2 = g^ v d^(j)d u 4>. Usual scalar fields are 
match to rj = 1 while rj = — 1 describe phantom fields. 
The homogeneous scalar field has the density p^ and 
pressure as 



P<p = \h 2 ^ 2 + v{4>) , 

P4> = ^H 2 cj)' 2 - V(<j>) 

and equation of state parameter obtain as 
1 - r 1 H 2 (/)' 2 /2V 



W0 = — = 



1 + riH 2 <i)' 2 /2V 



(49) 
(50) 

(51) 



Now the Friedmann equation take as (|Arianto et al 
20071) 



3/3 



(52) 



where (3 = f3 + -g^, . With this equation we obtain 
j>' = -2r,p(5£ + ?£) . (53) 



If this equation Substituting into the Friedmann equa- 
tio n, we can find th e potential of the scalar field 
as ( Arianto et al J 120070 



V = 3(3H 2 



H,4 
IT 



(54) 



where H = H((j>(t)). The equation of state find in 
following form 



-1 



-i-l P 



If 



A. 



(55) 



Relations (53) and (55) have main role in following con- 
siderations. We can see violation of the Lorentz invari- 
ance which has been described by existence of a vector 
field in the action, now it has combined in the dynam- 
ics of scalar field and equation of state by /3 parame- 
ter. This situation allows to consider phantom divided 
line crossing in LIV model. Now we should solve equa- 
tions, (53) and (55), to obtain dynamics of scalar field 
4> and the equation of state u)$. Therefore, we must 
define the Hubble parameter H(<p) and the vector field, 
In following, we choose a general case of the Hub- 
ble parameter H{<j>) and the vector field P{4>) and then 
consider crossing of phantom divided line, late time ac- 
celeration and Rip singularity scenario. 

3.1 Considering Late-Time Acceleration in LIV model 

For an acceleration universe, should be a > 0. We 
rewrite it to this form H'/H > — 1. For to consider 
this case, we must determine some equations. We study 
a general case for the vector field and the Hubble pa- 
rameter, note that these equations are a function of <fi 
(|Sadatian fc Nozari|[2008h 



H = H <f>' 



C 



71 > 2 



(56) 



Following, we just study a quintessence scalar field by 
rj = 1. By using relation (53), we calculate 



where 

A= \Ho(t -t )(-4:(m + 4:(mn + 2( 2 m 



(57) 



4 mn + 2 ran 2 ) + 4>o 



and by using relation (55), we can obtain 

uj^t) = -l + ^m<b n - 2 {t)(( + n) 2 . 

Following with relations (56) and (55,58), we find 
1 



m 2 < 



4(-l) n n - 2 (t)(C + 77) i 



n > 2, 



(58) 



(59) 



This equation describe a restriction in LIV parameters 
for a late time acceleration. 

Now with using equation (56) and (57) for obtain scale 
factor as 



i O™ + < 2) ( m(»+C)(»+C-2)-ff (t-t )+2«0 I = F -*0 



a(t) = a (t )e 



> + C-2) 



(n-2)m(n + C) 



9 



(60) 




in(2m(n + C)(rt + C-2)g (t-t ) + ^ 
e 1+C-2 



We should select a special space parameter, therefore 
we approach equation (60) by a Taylor series in special 
space parameter 

a(t) = .286504 + . 447663* + . 237821t 2 +.0352651t 3 (61) 
We can see in figure 5 acceleration evolution in late 




Fig. 5 Variation of scale factor a(t) for different values of 
t for n = 3 , ( = —2 and m = —0.1. The values of ( are 
determined by relation (59). 

time. Now, we study equation of state for determine 
crossing phantom divided line. By <p according to equa- 
tion (57), the equation of state obtain as 



where 



4 (C + n) s 

T H — m 

3 E 



(62) 



E= lH (t-t )(-4(m + 4(mn + 2( 2 m 
- 4mn + 2mn 2 ) + o ](~+c-2) ; 

which implicity has a dynamical behavior. This 
model lets us to choose a suitable parameter space to 
crossing phantom divided line. Moreover this parame- 
ter space must be compared by observational data. We 
insist in this model, a scalar field and a vector field 
together can explain crossing phantom divided barrier 
and late time acceleration. We can see in figure 6, the 




Fig. 6 ujff, vary for different values of the vector field by m 
and t for n = 3 and £ = — 2. Positive values of C, does not 
show a phantom divided line crossing. The values of ( are 
obtained with equation (59). 



crossing of phantom divided line for a dynamical equa- 
tion of state. Figure 6 perhaps explain why we are liv- 
ing in an era of tu < — 1. Note that f3(t) has the main 
role of Lorentz invariancc violation model. The dynam- 
ical equation for (3(t) has an interested mean: with a 
suitable fine tuning we can obtain a Lorentz violating 
cosmology correspond to observational data. 
Wc know many different models that explain phantom 
divided line crossing, but this model is special because 
it contains only a scalar field and a Lorentz violat 



ing vector field that checks the crossing (|Li et al.ll2008 



Libanov et al.l 120071 : iBertolami et al.1 120041 ) . However, 



two notice must be explained in this paper: First, in fig- 
ures 3 and 6 we can see, there are some sudden jerks of 
the equation of state. In some models that equation of 
state cro ssing the p hantom divided line, ui surge around 
— 1 (see (|Weill2006l ) and references therein). These jerks 
are certainly a signature of chaotic manner of equation 
of state. Second, in those figures we see, crossing of the 
phantom divided line appear at late-timc. This means, 
as second cosmological coincidence problem, requires 
extra fine-tuning in model parameters. 
For describing the physical reasons why the crossing 
of the phantom divide can occur in this model, we 
can state by a fitting coupling between a quintessence 
scalar field and other matter content may causes to 
a constant ratio of the energy densities of both con- 
texts which is reconcilable with an accelerated expan- 
sion of the universe or crossing of phantom dividing 
line. In this model, we have three kinds of energy- 
momentum: 1- standard matter, 2- Scalar Field as a 
nominee of Dark Energy and 3- energy-momentum con- 
tent rely on for Lorentz violating vector field. Here we 
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take standard matter has tiny contribution on the all 
energy-momentum content of the universe. For other 
two energy- momentum contents, it is possible to take 
the trigger mechanism to determined dynamical equa- 
tion of state. This means that we take scalar- vector- 
tensor theory containing Lorentz invariance violation 
which performs like the hybrid inflation models. In this 
condition, vector and scaler field play the roles of infla- 
ton and the waterfall field. Hence, it is reasonable to 
hope that one of them will eventually control to de- 
scribe inflation or accelerating phase and crossing of 
phantom divided line. 

3.2 Rip Singularity 

According to Rip singularity scenario that explained 
in previous section(page 5), we consider some types of 
Rip singularity in this LIV model. As shown in figures 
(5,6), we provide suitable conditions for considering Rip 
singularity scenario. Hence, correspond to equation of 
state, and using equations such as energy density, pres- 
sure and Hubble parameter, we can study some Rip 
singularity solutions. 

According our setup in this paper that selected ( = — 2 
and n = 3 in anzats (56) we have 



H = H 



f3 = mq 



Now with substituting above equations in energy den- 
sity and pressure relations, we can consider dynamical 
behaviors in these relations. Note that, we selected 
for ordinary matter p m = and p m — 0. Now after 
numerical calculations (figure 7), we obtain results as 
following: 

In finite time, scalae factor and energy density in- 
creases, but pressure decreases. Therefore type II sin- 
gularity (sudden) do not predict. But type III singu- 
larity appears because \p\ decrease so rapidly. In infi- 
nite time, Hubble rate tends to zero, therefore Little 
Rip singularity is impossible, but pseudo-rip singular- 
ity can occur. We emphasize this relations determine 
in original setup in this paper. This means by different 
selection in parameters of ( and to, all dynamic of rela- 
tions be change. Therefore, one can change setup and 
consider other options. However, in this situation, our 
work be limited. 



4 Comparison with Observational Data & 
Other Dark energy Models 



0.10 0.12 0.14 0.16 0.1K 



0.10 0.12 0.14 0.16 0.18 0.20 

t 

Fig. 7 Variation of pressure(up) and energy density(down) 
respect to time for a minimal coupling model(where £ = —2 
and n = 3) 



with precise observational data. Here we determined 
the data coming by SNe observations, the evolution of 
the Hubble parameter in the brane model. 

The modulus (i vs redshift z = a$/a—\ equation cor- 
responding to type la supernovae from the Supernova 
Cosmology Project ( Amanullah et al. 2010l: Uni 2013 ) 
is, 



H(z) = po + 5\gD L (z). 

The equation for the luminosity distance D^(z) as a 
function of the redshift in the FRW cosmology (FC) 
give as 



Df c = 



H 



(l + z) / h- l (z)dz, 



(63) 



I 1/2 



4.1 Analyzing SNe Data 



where h(z) = [O ro0 (l + zf + Ct D0 F(z) j 
Here, f2 m o is the total fraction of matter density, Quo 
the fraction of DE energy density, and Hq the current 
Hubble parameter. The constant value /zo conditional 
upon the chosen Hubble parameter: 



The parameters of the cosmological models may be con- 
sider from a explicit comparison of their predictions 



p = 42.384 - 5 log ft, ft = iVl00km/s/Mpc 
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The function F(z) = Pd(z)/pdo may be consider from 
the continuity equation 



PD - 3-g(p D ) = 0, 
a 



which can be reobtain as 
f PD(z) _dy_ 

J poo 9(y) 



-31n(l + z). 



(64) 



(65) 



For simplicity, we ignore the contribution of radiation. 
For cosmology on the brane (BC), Eq.(63) can be take 
as 

D? c = f h- 1 («)[l+<y/» 2 («)]- 1/3 (l+*) 1/a ds 

Ho Jo 

where for comfort the parameter 5 = Po/2 has been 
introduced. For the analysis of the SNe data one needs 
to obtain the parameter % 2 , which is given by 



2 _ V"^ [Pobs(Zi) — Pth(zi)Y 

xsn — 2^ — 



(66) 



where Uj is the corresponding la error. The parameter 
Po is free of the data points and, therefore, one has 
to do a uniform marginalization over p . Minimization 
with respect to po may be done by simply enlarging the 
Xsn with respect to po, 



Xsn = A — 2p B + plC 



(67) 



where 



[Pobs(Zt) - pth{Zi,Po = 0)]" 



Pobs{Zi) - Pth(Zi) 



The expression has a minimum for p = B/C at 

A - B 2 /C. 



Xsn 



One may minimize y.giy instead of X qm- 
with ( Nesseris fc Perivolaropoulosl 20051 ) and table 1, 
one obtains the 56.01% confidence level by A% 2 = 
X 2 — Xmin < 1-01 for the one-parameter or 1.9 for the 
two-parameter models, correspondingly, the 91.2% con- 
fidence level is determined by A% 2 = % 2 — Xmin < 3.78 
or 5.63 for the one- and two-parameter models. 



4.2 Dark energy Cosmology 

In this paper, we consider two models that explain 
dark energy in cont ent of the universe, but according 
( Bamba et al.l 120121) and references therein, there are 
many different models in this subject. Therefore, In 
this section for a worthy review, we briefly point out 
them. 

As we know some number of well-liked dark energy 
models, such as the ACDM model, Pseudo-Rip Lit- 
tle Rip and Little Rip scenarios, the phantom and 
quintessence cosmology with the four types of the sin- 
gularities and non-singular universes satisfied with dark 
energy have been considered. 

Investigations have been shown the ACDM model 
and various cosmological observations to endure the 
bounds on the late-time acceleration of the universe. 
In this regards, researchers have considered a fluid de- 
piction of the universe in which the dark fluid has a 
popular form of the Equation of state with the inho- 
mogeneous. They have shown that all the dark energy 
cosmology can be describe by some different fluids, also 
determined their properties. It has also been display 
that at the present stage the cosmological evolutions of 
all the dark energy universes can be like to that of the 
ACDM scenario, and hence those models are satisfied 
with the cosmological observations. In other hand, they 
have studied the equality of appreciate of different dark 
energy models, this means, single and multiple scalar 
field theories could be portray. 

However, we emphasize the role of cosmography in this 
review. It is a essential tool because it lets, in doctrine, 
to separate among models without a seizures but lay- 
ing on restriction coming from data. However, the goal 
criticism to this approach is depended to the expansion 
of the Hubble parameter. In mostly, observations could 
not be enlarged at any red-shift and, in some of occa- 
sions, are not enough to track models up to early eras. 
Finally, the future observational data campaigns must 
fix the condition elimination the degeneration incipient 
at low red-shifts and let a well insight of models. 



z 


H bs(z) 

km s _1 Mpc -1 


OH 




km s" 1 Mpc -1 


0.170 


83 


8 


1.530 


140 


14 


1.750 


202 


40 



Table 1 Hubble parameter in contrast to redshift data. 
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5 Summary 



Appendix 1: Proof of stability 



Several evidences from observational data have shown 
that our universe is in the accelerated expansion era. 
In this paper, we have shown this reality in two dif- 
ferent models: A non-minimally coupling scalar field 
identified on the brane and a scalar field coupling min- 
imally to gravity in LIV model. According complex of 
dynamical equations, we have limited our research to 
some special form of non-minimal coupling and scalar 
field potential. Then we have studied special form of 
time evolution for the scale factor and a scalar field. In 
brane world model, we have determined a equation that 
explain situations where an accelerating expansion of 
the universe is possible. As a consequence, we have ob- 
tained which non-minimally coupling scalar field iden- 
tified on the brane is a suitable candidate for late-time 
expansion of universe. In LIV model, we have consid- 
ered a new framework for crossing of phantom divided 
line with equation of state, that combined a violation of 
Lorentz invariance in a cosmological model. However, 
we have explained with a suitable choice of parameter 
space, it is possible to have phantom divided line cross- 
ing just with a Lorentz in variance violating vector field 
and a single scalar field (INesseris fc Perivolaropoulos 
2007t IVikmanl 120051) . In this regard, existence of a 
Lorentz invariance violating vector field prepare a setup 
for crossing phantom divided line just with a minimally 
coupling scalar field. In other view point, there is the 
possibility of a Rip singularity by suitable tuning in the 
parameters of models. 

Comparison of non-minimal model with WMAP data 
takes more precise restricts on the values of non- 
minimal coupling. A detailed comparison between 
this results and some results of previous considerations 
shows that our constraints for non-minimal coupling 
with exponential potential are steadfast by holographic 
dark e nergy and also with result of warp ed DGP brane 
model (jNesseris fc Perivolaropoulosll2007l ). 



Now we consider the stability of our model. In or- 
der to consider the classical stability of our model, we 
determined the manner of the model in the uj — uj 1 plane 
where u' is the d erivative of to with th e logarithm of the 
scale factor (see (|Arianto et al.l 120071 ) for a comparable 
analysis) 



du> 



duj dt 



d In a dt d In i 



H' 



(68) 



The sound speed declares the phase velocity of the in- 
homogeneous perturbations of the field. Hence, we take 
the function c a as 



c 2 = l 

( Y7 . • 



(69) 



If the matter is determined as a perfect fluid, this equa- 
tion could be the adiabatic sound speed of this fluid. We 
emphasize with scalar fields that do not comply perfect 
fluid form necessarily, this parameter is not actually a 
sound speed. The conservation of the total (scaler and 
vector ) field given with (45). As we view, that relation 
implicity depended on /? and H form. Therefore, we 
can take anzats (56) with substituting in (45) as 



dptotal 

dt 



+ HPtotal +P total) = 



(70) 



Since the dust matter complies the continuity equation 
and the Bianchi identity keeps legally , total energy 
density fulfill the continuity equation. From that equa- 
tion, we have 



Pde 



-3/Ode(l + w de 



(71) 



where de means dark energy. By using equation of state 
Pde = UdePde, we have 



Pde = upde + U de p de 

Hence, the function c a 2 can write as 

•2 ^de 



-3(1 + LJde) 



(72) 



(73) 



In this condition, with suitable choice of £ and n in (56), 
which we take in pagc(10), we calculated the lo — u' 
plane is divided into four regions defined as 



U!de > — 1, > 3w(l + Ul) 



Ca' > 



II : uJde > -1, w' < 3w(l + us) 



< 



III : code < -1, w' > 3w(l + Ul) Ca < 
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IV: ui de <-l, ui' <3oj(1 + cj) => c a 2 > (74) 

As one can see from these equation, the regions I and 
IV have the classical stability in this model. 



Appendix 2: initial conditions for numerical calcu- 
lations in both models 

A : In the brane world model 
In the induced gravity model we should add the extra 
condition 

p> \6a{H 2 + — )| 

Hence we need both high energy and weak coupling. 
In other hand, the last condition and the continuity 
equation hints that P = — p for H ^ 0. However that 
unlike 4D general relativity, we do not necessarily claim 
p = constant. In other view point, initial values for 
parameters £ and ft given from equation 



1 (3^ - l) 2 
\Q2v{2v- 1) 



and 



iv- 1 



± 



{Zv — l) 2 



1§&(2v- 1) 



B : In Lorentz invariance Violation Model 
The dynamical attractor of the cosmological system 
has been occupied to make the late time manners of 
the model uncaring to the initial condition of the field 
and therefore alleviates the fine tuning problem. In 
quintessence models, the dynamical system has track- 
ing attractor that makes the quintessence develops by 
tracking the equation of state of the background cosmo- 
logical fluid so as to relieving the fine tuning problem. 

Hence, in order to make viable Lorentz violation 
model, we need that the coupling function /3 and the 
potential function V must satisfy the condition 



and 



VV^ 1 Mi 
VI + 2V,/V > 



respectively. 



li 
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